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Abstract. A radio continuum survey at medium Galac- 
tic latitudes with the Effelsberg 100-m telescope is being 
carried out at a centre frequency of 1.4 GHz in total power 
and linear polarization. Areas up to ±20° of Galactic lati- 
tude are now being observed at a sensitivity of 15 mK Tb 
in total intensity and 8 mK Tb in linear polarization with 
an angular resolution of 9f 35. This paper describes the ob- 
serving and reduction technique applied which results in 
absolutely calibrated maps. The methods are illustrated 
by examples of images from the survey. 
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1. Introduction 

The interaction between the Galactic magnetic field and 
interstellar clouds leads to a variety of radio emitting 
structures, shells, filaments, and loops. Many of these fea- 
tures are linearly polarized. They extend up to high Galac- 
tic latitudes but they are very faint, an order of magni- 
tude fainter than the large-scale diffuse Galactic emission. 
However, the study of these structures may provide im- 
portant clues for the understanding of the "disk-halo con- 
nection" . The interpretation favours chimneys, fountains, 
Parker loops or superbubbles created by OB associations 
to transfer material from the disk into the halo. They will 
deform and compress the Galactic magnetic field. If these 
structures account for the rather smooth nonthermal back- 
ground, they have to be sufficiently numerous and large. 

The Galactic region above b = +4° and below b = —4°, 
where the just mentioned structures play an important 
role, has never been studied in a systematic way. Large- 
scale surveys covering this area either suffer from low an- 
gular resolution or they miss the diffuse emission as it is 
the case for the 1.4 GHz VLA-survey (Condon et al. |1996Q . 
A combination of 1.4 GHz data from the Effelsberg tele- 
scope with high angular resolution data from the VLA is 



of particular importance to separate compact background 
sources from faint extended Galactic structures. An ex- 
ample has already been shown and discussed by Fiirst et 
al. ( |1998D . 

One of the milestones in the research of Galactic mag- 
netic fields is due to Brouw & Spoelstra ( |1976| ) and Spoel- 
stra (1984). They have surveyed the northern sky in lin- 



ear polarization with the Dwingeloo 25-m telescope at 
1.411 GHz at an angular resolution of about half a degree. 
They also compiled maps of the linear polarization at fre- 
quencies bet ween 4 08 MHz and 1.411 GHz. At 2.7 GHz 
Junkes et al. ( 1987 ) published a survey of the linear polar- 
ization of the Galactic plane (|6| < 1?5 and 4? 9 < I < 76°). 
Recently, the southern Galactic pla ne Q6 | < 5°) was sur- 
veyed at 2.4 GHz by Duncan et al. ( 1995 ). Aside from the 
Dwingeloo data (see the discussion in Sect. 5), all these 
polarization surveys are constrained to the Galactic plane. 
On th e low frequency side, at 327 MHz, Wieringa et al. 
(1993) observed small-scale variations in polarization at 
high latitudes with the Westerbork synthesis telescope. 

The knowledge of the magnetic field in the Galactic 
halo even in the immediate vicinity of the Galactic plane 
is limited, because of the lack of observations of the linear 
polarization at the proper angular resolution and sensitiv- 
ity. To fill this gap new observations at 1.4 GHz are being 
carried out to cover the entire Galactic plane at medium 
Galactic latitudes (|6| < 20°) visible at Effelsberg. This 
survey will reach a sensitivity at total intensity close to 
the confusion limit (~ 15 mk Tb) and an even higher sen- 
sitivity (~ 8 mK Tb) at linear polarization. 

We discuss the general survey parameters and observ- 
ing method in Sect. 2. An absolute calibration method 
for the total intensity data is described in Sect. 3. The 
method developed to overcome the usual problem of po- 
larimetric observations, the instrumental polarization, is 
demonstrated in Sect. 4. In Sect. 5 we propose a procedure 
to adjust the polarization data to an absolute temperature 
scale. 
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2. Observing method and data reduction 



Table 2. Major calibration sources 



The observations of the survey are carried out with the 
two-channel 1.3-1.7 GHz receiver installed in the primary 
focus of the Effelsberg 100-m telescope. The receiver is 
equipped with cooled HEMT amplifiers recording the left- 
and right-hand circularly polarized component simultane- 
ously (see Schmidt & Zinz 1994). These amplifiers are of 
very high stability, and total intensity observations are 
limited by confusion which is reached after about 2 s of 
integration time. An IF-polarimetcr converts the circu- 
larly polarized component into linear components (Stokes 
parameter U and Q). Stokes parameter / is obtained by 
simply adding both circularly polarized components. 

The method of observation is to scan each field along 
constant Galactic latitude as well as along constant Galac- 
tic longitude. Each field was measured at least twice. The 
telescope parameters relevant to the 1.4 GHz survey are 
listed in Table 0. 



Source 


I 


b 


S^lcm/Jy 


% Pol. 


PA 


3C286 


56.5 


+80.7 


14.4 


9.3 


32 


3C138 


187.4 


-11.3 


9.5 


7.9 


169.9 


3C 48 


134.0 


-28.7 


15.9 


0.4 


139.7 



Table 1. Parameters for the survey 



Parameter 



Value 



Centre Frequency 
Bandwidth 

Integration time per point 
System temperature 
Typical rms noise in total power 
Typical rms noise in polarization 
Tb / S 

Telescope beamwidth 
Typical scan length 
Scan interval 
Scanning velocity 



1.4 GHz* 
20 MHz* 
2 s 
26 K 
15 mK 
8 mK 

2.12 ±0.02 K/Jy 
9 .'35 ± 0.'04 
10° 
4' 
4'/s 



* sometimes interference makes a shift of the centre frequency 
and/or a reduction of the bandwidth necessary. 



Daytime observations cause inacceptable distortions in 
both the total intensity and the polarization data due to 
contribution of solar emission to the far sidelobes of the 



telescope (Kalberla et al. 1980). Therefore all observations 
are carried out at nighttime, when additional ionospheric 
Faraday rotation effects are minimal and included by cor- 
rections inferred by observations of calibration sources. 
Calibration sources were 3C286, 3C138 and 3C48, where 
3C286 served as the primary calibrator both in total in- 
tensity and polarization. The scale accuracy for total and 
polarized intensities is better than 5%. As seen from their 
coordinates, calibration sources have large angular separa- 
tions on the sky, so that at least one of them is accessible 
at all times. Their polarization properties are listed in Ta- 
ble]^ and are taken from the list of Tabara & Inoue ( 1980| ). 

In Fig. [I] we show a beam map at 1.4 GHz observed 
in the Az/El-system on 3C123, which at that frequency 
shows no measurable variation with elevation. Compact 



Fig. 1. 1.4 GHz beam map covering an area of 1?1 x 1?1. 
Contours run in steps of 3 dB. The lowest contour is at 
—30 dB and the highest contour is at —3 dB of the maxi- 
mum. 

background sources in the field of 3C123 have been re- 
moved. The dynamic range exceeds 30 dB. The maximum 
of the first sidelobes, which are enhanced by the four sub- 
reflector support legs, are at a level of about —20 dB. 
The cross-polarization lobes of Stokes parameters U and 
Q are dominated by a main beam component, which varies 
with parallactic angle. A correction of that instrumental 
response is described below. 

The first stage of the data reduction made use of the 
standard "Toolbox" procedure (CONT2, von Kap-herr 



1977) for continuum and polarization observations with 
the Effelsberg telescope. From the tabulated scans cali- 
brated maps were computed for all the four data channels 
(two times /, U and Q) recorded. A linear baseline us- 
ing data points at the end of each scan was subtracted 
from each channel. All later stages of the data reduc- 
tion are based on the NOD2 program package (Haslam 



1974). Spiky data were removed. Scanning effects were 



suppressed by using the method of unsharp masking de- 



veloped by Sofuc & Reich (1979). Finally, the two maps 



observed in orthogonal directions for each of the Stokes 
parameters have been added using the PLAIT program 



(Emerson & Grave 1988), which reduces "scanning ef- 
fects" from individual distortions or by the baseline setting 
procedure in the case of an emission structure located at 
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the edge of a map by appropriate weighted addition of the 
Fourier transforms of both maps. The polarization U and 
Q maps were corrected for instrumental polarization as 
described in Sect. 4 before combining the maps observed 
in orthogonal directions. 



3. Absolute calibration of the Effelsberg total 
intensity maps 

From the observation and the method of data reduction it 
is obvious that neither the total intensity nor the polariza- 
tion maps are on an absolute temperature scale. The total 
intensity maps are calibrated to an absolute scale using the 
1.4 GHz northern sky survey by Reich (1982) and Reich & 
Reich (1986) carried out with the Stockert 25-m telescope. 
This procedure has already been described by Reich et al. 



(199C) when calibrating the Effelsberg 1.4 GHz Galactic 
plane survey. The Stockert survey has a scale accuracy of 
5% and an uncertainty of the absolute zero-level of 0.5 K. 

Briefly, the Effelsberg data and the Stockert data are 
convolved to an angular resolution slightly exceeding that 
of the Stockert 1.4 GHz survey (HPBW 36') and the dif- 
ference between the two data sets is added to the Effels- 
berg data. This procedure also improves the original Ef- 
felsberg maps since large-scale distortions by atmospheric 
and ground radiation variations are removed. 

However, still existing faint baseline effects in the low 
resolution Stockert data are also added to the higher reso- 
lution Effelsberg data. Such distortions were found to exist 
on small scales but not on the large-scale background com- 
ponent, which is the missing component of the Effelsberg 
maps. A modified version to calibrate the Effelsberg maps 
absolutely was introduced: 

• Convolve the Effelsberg map to the resolution of the 
Stockert map (36'). Decompose this map into a large- 
scale component (Eff.back) and a map with small-scale 
structures using the "background filtering method" 
(Sofue & Reich 1979 ) with a 3° Gaussian beam for 
smoothing. 

• Subtract the 2.8 K isotropic background component 
from the Stockert data and convert the temperature 
scale from full beam into main beam (Xmb/7fb = 
1.55). 

• Decompose the Stockert data to obtain the corre- 
sponding large-scale component (Sto.back) by using 
the same filtering parameters as for the Effelsberg 
data. 

• Adjust Eff.back to Sto.back as described above, 
namely, add the difference of both large-scale compo- 
nent maps to the Effelsberg large-scale data and add 
the 2.8 K isotropic background component. The result- 
ing map is ES.back. 

• The final Effelsberg map (Eff.final) is given by: 

Eff.fmal = Eff.orig + ES.back - Eff.back. 



Fig. 2. A small section from the survey, towards £ ~ 50°, 
illustrates the calibration of the total intensity data to the 
absolute temperature scale. The panels, from top to bot- 
tom, display the Effelsberg (contours start at 100 mK and 
are plotted in 150 mK steps) and Stockert 1.4 GHz (con- 
tours start at 4500 mK in 100 mK steps, full beam bright- 
ness scale) measurements and the combination of the two, 
respectively. Contours for the combined map (bottom) run 
from 5500 mK in steps of 150 mK. 
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We use the modified method only when there are non- 
negligible scanning effects in the source component of the 
Stockcrt data. A sample region (towards I ~ 50° ) is given 
in Fig. U to illustrate the effect of the absolute calibration. 

4. Instrumental Polarization 

One major goal of the new 1.4 GHz survey is the sen- 
sitive mapping of polarized emission out of the Galactic 
plane. However, a problem is given by the relatively high 
instrumental polarization mainly introduced by the cooled 
broad-band polarization transducer and hybrid of the L- 
band Effelsberg receiver. A change of the bandwidth and 
centre frequency has a significant effect on the instrumen- 
tal U and Q components. Other instrumental effects due 
to the antenna and feed characteristics, residual ellipticity 
of the polarimeter response and variations with time add, 
but are too small to be separated from the main effect. 

An attempt was made to minimize these instrumental 
effects to a residual effect of the order of 1%. We assume 
that the instrumental components U- ms t and Q ms t scale 
with the total intensity /. When observing in an astro- 
nomical coordinate system the instrumental components 
depend on the parallactic angle (f>. A procedure was de- 
veloped to correct for this effect. From the observations 
of the polarized calibration sources 3C286 and 3C138 at 
different parallactic angles, correction factors /tj and /q 
are calculated in such way that the nominal percentage 
polarization and the polarization angle are obtained. 

The algorithm of the elimination of the instrumental 
effects is as follows: Suppose t/ ac t and Q ac t are the in- 
trinsic values of a source in an astronomical coordinate 
system. U oos and Qobs are the observed values depending 
on the parallactic angle <f>. The instrumental values C/i ns t 
and Qi ns t are: 



U iast (4>) = U bs(4>) - U act 

Qinst(</>) = Qobs ((/>)— Qact 



(1) 



The instrumental polarization angle PAi nst and polar- 
ization intensity -Pimst are: 



PA inst {(l>) 



1 



■ tan 



£Anst(0) 
<9inst(^) 



(2) 



The correction factors /tj and /q are then: 



M<P) 
/qWO = 



PI instil) 

lobs 
P/instM 

^obs 



sin [2(PA iast (4>) - 0)] 
cos [2(PA inst ((j>) - (j>)] 



(3) 



In Fig. [|we show, as an example, the data for /tj and 
/q from observations of 3C286. /tj and /q are determined 
for each session. It is obvious that /tj and /q vary with 
parallactic angle, which needs to be taken into account 



when correcting a survey map. According to the intro- 
duced procedure each pixel / of a map is multiplied with 
the appropriate /tj and /q. The resulting U and Q com- 
ponents of the instrumental polarization are transformed 
into parallactic angle corrected components for each pixel 
of an U and Q map and subtracted from the observed U 
and Q. The procedure corrects for most of the instrumen- 
tal effects. The corrected U and Q maps are then used to 
calculate the polarized intensity PI and the polarization 
angle PA. 



Fig. 3. Variation of the calculated /tj and /q factors with 
respect to the paralactic angle for 3C286 during one ob- 
serving session. Error bars are the standard deviations of 
the factors. 

It has been found that the correction effect on the 
large-scale polarization emission is not significant. Strong 
sources, however, cause distortions which are clearly min- 
imized by the described procedure. In Fig. [| we illustrate 
the effect of our correction procedure. We note that in- 
strumental effects do not always act in a way to increase 
the observed polarization. Depending on the arrangement 
and variations in U and Q it is equally possible to observe 
the reverse of this effect. 

5. Procedure for the absolute polarization 
adjustment 

As was mentioned earlier the large-scale information in the 
polarization maps is lost due to the baseline fitting pro- 
cedure. In addition, residual ground radiation effects may 
show up in view of the large sizes of the observed maps. As 
in the case of total power maps, the polarization emission 
needs to be adjusted to an absolute level. However, there 
exists no complete and regularly gridded low resolution 
polarization survey which is absolutely calibrated. Fortu- 
nately, a large set of linear polarization observations at 
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1.4 GHz obtained with the Dwingeloo 25-m telescope ex- 
ists (Brouw & Spoelstra 1976). These data are absolutely 
calibrated and corrected for all kinds of instrumental ef- 
fects and therefore match the requirements for this task. 
However, the data are not on a regular grid and are signif- 
icantly undersampled. This data set has to be regridded 
to be used to adjust our polarization measurements. The 
polarized intensity and polarization angle data from the 
Dwingeloo survey have been provided by Dr. Spoelstra in 
digitized form. 

Using the PI and PA values from this measurements 
we calculated U and Q values and regridded the under- 
sampled Dwingeloo data on the grid of the Effelsberg 
maps. However, in some regions data points are too sep- 
arated (2° or more) for an absolute calibration, and in 
other weakly polarized regions an S/N-ratio of 2 or less 
does not allow a proper adjustment. 

We tried two methods to interpolate the Dwingeloo 
data on the same grid as the Effelsberg maps. One way 
to do this is an interpolation between the data points us- 
ing a "cubic-spline interpolation" . However, we found that 
a cubic-spline interpolation introduces distortions at the 
corners of the maps and data from a much larger area 
must be used to avoid this problem. Moreover, a single 
high intensity data point affects surrounding low intensity 
data up to a large distance. Data points must be weighted 
with respect to the distance. Hence, we found that a cubic- 
spline interpolation is inadequate for most of our regions. 

A successful method is to weigh data points by their 
distances to the required map element. We used the ap- 
proach exp(— a TZ) in which a is a constant. We calculate 
for each pixel of an Effelsberg map the Dwingeloo U and 
Q data within a radius, TZ, contributing with a weight as 
given above. For our case we found a value of 6° for TZ 
with a = 1 to be satisfactory. 

The reconstructed U , Q and PI Dwingeloo maps of a 
test region at the 4' grid of the Effelsberg maps are given 
in Fig. |. 

Figure ^| shows the U and Q maps of the original Ef- 
felsberg measurements of an area, which is a small section 
of the Dwingeloo map shown in Fig. |^. In the higher reso- 
lution Effelsberg maps numerous small-scale polarization 
structures are visible, which are smoothed out by the large 
Dwingeloo beam. 

The two data sets are combined as follows: The Ef- 
felsberg map is convolved to the Dwingeloo beam (36') 
and subtracted from the Dwingeloo map. The difference 
is added to the original Effelsberg map. Figure and Fig. | 
show the combination for U and Q and the correspond- 
ing PI map. In these figures the small-scale structures are 
much less pronounced due to the addition of the strong- 
large-scale polarized emission which varies in the range 
from 200 mK to 800 mK across the map (Fig. ||). 

The absolutely calibrated U and Q maps may be de- 
composed into small-scale and large-scale features by stan- 
dard methods for total intensity maps and relative PI 



maps can be calculated. Some examples have been given 



by Uyamker fll997| ). 

Plots of the polarization vectors demonstrate the ef- 
fect of adjusting the Effelsberg data to an absolute level. 
The polarization angle maps are presented as vector plots 
in Fig. ^. Numerous small-scale structures are visible in 
the original Effelsberg map. The polarization angle varies 
largely across the maps. However, the electric field vectors 
are almost constant in the combined Effelsberg-Dwingeloo 
map. The data of this map are the same as those for Fig. |[ 
but grey-scale representation is much more sensitive to 
small variations. 

6. Final remarks 

We described the observation and reduction technique for 
a sensitive 1.4 GHz continuum and polarization survey 
with the Effelsberg 100-m telescope. We illustrated by an 
example the absolute calibration of the total intensity Ef- 
felsberg data by making use of the 1.4 GHz Stockert survey 
and showed the effect on the Effelsberg polarization maps 
after adjusting the large-scale polarized emission with the 
Dwingeloo 1.4 GHz measurements. 

The methods used and introduced in this work have 
been found to be appropriate to conduct a 1.4 GHz survey 
which is absolutely calibrated both in total and polarized 
intensity and corrected for instrumental polarization down 
to a level of 1%. 

Acknowledgements. We are very grateful to Dr. Titus Spoel- 
stra for providing the Dwingeloo polarization data in digital 
form. We thank Dr. W.N. Brouw for helpful comments on the 
manuscript. 
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Fig. 4. A sample region to illustrate the effect of the procedure to eliminate the instrumental polarization. The upper 
panel displays the same region before (left) and after (right) the correction. Contours show the total intensity starting 
from —50 mK with 50 mK intervals. The strong source around I ~ 47° is 3C386 with a flux density of 6.5 Jy. Its 
percentage polarization is 2.7% after the procedure is applied. In the lower panel the polarized intensity data of the 
same regions are displayed in grey scale. Contour levels run starting from 10 mK with 30 mK steps. Again the original 
data are at left and corrected data are at right. The plotted electric field vectors are scaled to the polarized intensity 
and astronomical coordinates such that 100 mK correspond to a vector of length 1.'5. 
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Fig. 5. Polarization maps reconstructed from the Dwingeloo data as explained in the text. The panels, from top to 
bottom, display the Stokes U and Q maps and polarization intensity, respectively. The area marked with dashed lines 
is the region observed with the Effelsberg telescope and used to demonstrate the method of absolute calibration for 
the polarization data. 
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Fig. 6. Polarization U and Q maps from the Effelsberg 1.4 GHz survey. Contours are plotted starting from the lowest 
value of the wedge and in steps of 30 mK Tb- Contours starting from zero are plotted in white with the same steps. 
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Fig. 7. Polarization U and Q maps after calibrating the Effelsberg polarization maps to absolute temperature scale. 
Contours are plotted starting from the lowest value of the wedge and in steps of 30 mK Tb . For U map the contours 
starting at 200 mK Tb and for the Q map the contours starting at 400 mK Tb are plotted in white. Contour steps 
are always 30 mK Tb . 
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Fig. 8. Polarization intensity maps before (top) and after (bottom) calibrating the Effclsbcrg polarization maps to 
absolute temperature scale. Contours are plotted starting from the lowest value of the wedge and in steps of 30 mK 
Tb- For the upper panel contours starting at 120 mK Tb and for the lower panel contours starting at 450 mK Tb are 
plotted in white. 



12 



B. Uyamker et al.: A 1.4 GHz radio continuum and polarization survey at medium Galactic latitudes 



Fig. 9. Electric field vectors before (top) and after (bottom) calibrating the Effelsberg polarization maps to absolute 
temperature scale. The lower panel is essentially the same as Fig. 8. Due to the high polarization intensity of background 
polarization small-scale variations from the Effelsberg observations are hidden. Electric field vectors are scaled to the 
polarized intensity such that 100 mK correspond to a vector of length 3' for the upper panel and 1' for the lower panel. 
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